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Motivated by recent progress in development of cryogenic memory compatible with 
single flux quantum (SFQ) circuits we have performed a theoretical study of magnetic 
SIsFS Josephson junctions, where 'S' is a bulk superconductor, V is a thin super- 
conducting film, 'F' is a metallic ferromagnet and T is an insulator. We calculate 
the Josephson current as a function of s and F layers thickness, temperature and 
exchange energy of F film. We outline several modes of operation of these junctions 
and demonstrate their unique ability to have large product of a critical current Iq and 
a normal-state resistance Rn in the ir state, comparable to that in SIS tunnel junc- 
tions commonly used in SFQ circuits. We develop a model describing switching of the 
Josephson critical current in these devices by external magnetic field. The results are 
in good agreement with the experimental data for Nb-Al/AlO^-Nb-Pdo.ggFeo.oi-Nb 
junctions. 

PACS numbers: 74.45,+c, 74.50.+r, 74.78.Fk, 85.25. Cp 
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Practical applications of superconducting digital circuits were significantly limited by the 
relatively low capacity of superconducting memory. This motivated initial proposals to use 
superconductor /ferromagnet (S/F) hybrid structures as basis for the development in cryo- 
genic magnetic Random Access Memories (RAMs) 1-3 . Following the first experimental real- 
ization of SFS Josephson junctions 4,5 , much attention was paid to realize Josephson devices 
with complex magnetic barriers allowing switching between high and low critical currents. 
A number of different device structures were considered^— based either on superconducting 
spintronics effects or singlet-triplet switching within the magnetic barrier. However, these 
approaches were based on structures with reduced characteristic voltage I c Rn of junctions. 

Recently, successful realization of switchable Nb-Al/AlO^-Nb-Pdo.ggFeo.oi-Nb junctions 
was reported . These junctions are of SIsFS type, i.e., a serial connection of the Sis 

tunnel junction and sFS sandwich. SIsFS structure has high characteristic voltage, I c Rn, 
due to the presence of tunnel barrier 'I'. At the same time the whole structure behaves 
as a single junction with respect to an external magnetic field H ext and magnetic flux $ 
penetrating into the structure, since the intermediate layer s is too thin to screen magnetic 
field. As a result, the magnetic field entering the Pdo.99Feo.01 layer will modify its effective 
magnetization, facilitating the critical current control of the whole double-barrier SIsFS 
structure. According to Ref.— , effective magnetization in the dilute Pdo.99Feo.01 is controlled 
by Fe-rich PdsFe nanoclusters, which can be easily reordered by a weak magnetic field. 
The purpose of this paper is to develop a theory describing various modes of operation of 




FIG. 1. Schematic design of SIsFS Josephson junction. Solid line demonstrates typical distribution 
of pair potential A along the structure. It reaches bulk values in the S electrodes, is suppressed 
in the middle s layer and vanishes in the ferromagnetic region F. The characteristic length scales 
are also marked in the figure: Az, is the London penetration depth and £s is the coherence length 
typical for niobium. 
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SIsFS Josephson devices. We propose a quasiclassical model for a double-barrier Josephson 
structure in which magnetic and spin states of the F-film in an sFS part of the structure 
can change the properties and even the ground state of its Sis part. As a result, an external 
magnetic field may switch the junction from a superconducting to a resistive state or may 
transform a conventional current-phase relation to an inverted one. We compare the results 
with experimental data recently obtained for for Nb-Al/AlOz-Nb-Pdo.ggFeo.oi-Nb junctions. 

We consider the multilayered structure presented in Fig.l. It consists of two supercon- 
ducting electrodes 'S' separated by a tunnel barrier T, an intermediate thin superconducting 
film V and a ferromagnetic layer 'F\ To describe the supercurrent transport in the struc- 
ture, we assume that the conditions of a dirty limit are fulfilled for all metals. We also 
assume that all superconducting films in the structure are made from identical materials, 
i.e., they can be described by the same critical temperature, Tc, and coherence length, 
£s — (Ds/2nTc) 1 ^ 2 , where D$ is the electronic diffusion coefficient. The tunnel barrier I 
and the sF and FS interfaces are characterized, respectively, by the following parameters 
7b/ = RjA/^sPs, Ibfs = Rfs-A/^fPf, and 7 = Ps£s/£fPf- Here A, Ri and R FS are the 
area and the resistances of the interfaces, £$,f and ps,F are the coherence lengths and normal 
state resistivities of S and F materials, respectively. Under the above conditions the Joseph- 
son effect in the SIsFS junctions can be described by solving the Usadel equations^"— with 
Kupriyanov-Lukichev (KL) boundary conditions^ at Is, sF and FS interfaces and with the 
bulk pair potential in the depth of S-electrodes. 

The formulated above boundary problem has been solved numerically. The results are 
summarized in Figs. 2-4, where various modes of operation of the structure are defined ac- 
cording to chosen materials and layer thicknesses. These modes are clearly defined by the 
dependencies of characteristic voltage IqRn on thickness of intermediate superconductor, 
L s , and ferromagnetic, L F , layers (see Fig 12]- FigfjJ. 

Mode (1) in Fig. 2. If the thickness of the middle s-electrode, L s , is much larger than 
the critical thickness of s layer, L s c, which separates the different modes of operation in 
SIsFS structures, the pair potential A in the s layer is close to that of bulk material. Note 
that the critical thickness L s q in an sN (sF) bilayer at a given temperature is generally 
defined as a minimal thickness of an s-layer when superconductivity still exists. In the mode 
(1) the structure can be considered as a pair of Sis and sFS junctions connected in series. 
Therefore, the properties of the structure in parameter range (1) are almost independent 
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FIG. 2. Characteristic voltage IcRn of the SIsFS structures versus thickness of the F-layer Lp for 
different thicknesses of the middle superconducting film L s at T = 0.5Tc- Short-dashed straight 
line shows the IcRn product of the tunnel SIS junction at the same temperature. Interface 
parameters: jbi = 1000, jbfs = 0.3 and 7 = 1 at the sF and FS interfaces. 

on the thickness L s and are determined by the junction with smallest critical current. It is 
seen from Figj2]that for the given parameter set, T = 0.5T C , H = 10T C , 7 = 1, 7b/ = 1000, 
Ibfs — 0.3, the critical thickness of the s layer, L s c, is close to 3£s- 

Mode ( la) in Fig\^ In the ordinary case of Icsis IcsFS, the behavior of the structure 
coincides with that of conventional SIS junction with one important distinction - the sFS 
junction can turn the SIsFS structure into a 7r-state. At the same time, other properties 
like high IcRn product and sinusoidal current-phase relation are preserved in the 7r-state. 
Therefore the structure can be called switchable 0-n SIS junction. 

Mode (lb) in Fig\^ Another limiting case is realized for large Lp values and high ex- 
change fields H. Namely, the structure transforms into a standard SFS-j unction without 
any influence of tunnel barrier. 

Mode (2) in FiglE The absence of superconductivity in the s-electrode in the opposite 
case (L s L s c) leads to formation of the complex -InF- weak link area, where n marks 
the intermediate s film in the normal state. It results in much smaller critical current value 
Ic, with the magnitude close to that in well-known SIFS junctions^. The dependence of 
Ic on the thickness L s is weak due to large decay length in the n-region with suppressed 
superconductivity. 
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FIG. 3. The temperature dependence of characteristic voltage IcRn of SIsFS structure for dif- 
ferent values of exchange field H in the F-layer. The short-dashed line demonstrates the depen- 
dence characteristic for a conventional tunnel SIS Junction. It is seen how the exchange field 
H shifts the effective critical temperature Tq, corresponding to the switching of the s-layer from 
the superconducting state to the normal one. The circles show IqRn measured in Nb-Al/A10 x - 
Nb-Pdo.ggFeo.oi-Nb junctions^, proving the existence of effective critical temperature Tq in these 
samples. 

Mode (3) in Fig\^ Conversely, in the intermediate case (L s « L s c) the properties of 
the structure are extra sensitive to variations of the decay lengths parameters. Within the 
considered intermediate thickness range the system may transform from the mode (1) to 
the mode (2). Moreover, in this situation the system is sensitive to the F-layer parameters 
(thickness Lp and exchange field H), since these parameters control the suppression of 
superconductivity in the sF bilayer. 

This sensitivity allows one to change an operation mode by changing the parameters 
such as effective exchange field H and temperature T. Note, that depending on the domain 
structure of a ferromagnet and morphology of the F-film it might be possible to control the 
effective exchange field H. 

FigfS] demonstrates the temperature dependence of the critical current in the structures 
with thickness around critical one ( L s c = 3£s) for different values of exchange field H. 
These structures are characterized by the existence of the effective critical temperature T£ 
which corresponds to the appearance of superconductivity in middle s-layer and, correspond- 
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ingly, to an exponential growth of the current. Therefore, T c may significantly shift during 
remagnetization of the system (due to changing of H, as pointed out above). Thus, system 
can exist either in the superconducting or in the normal state depending on the history 
of the application of a magnetic field. On the other hand, from point of view of practical 
applications, the 0-tt SIS mode (la) seems more relevant. 




FIG. 4. The dependence of characteristic voltage IcRn on F-layer thickness Lp in the SIsFS 
structure with the s-layer in the superconducting state. Inset shows the current-phase relation in 
the vicinity of the first — n transition. Switching from to n state in the mode (la) preserves 
the value of critical current Iq as well as characteristic voltage IcRn- 

FigJH demonstrates that change of F-layer thickness Lp leads to 0-7r transition. The 
system can be switched into a 7r-state keeping the value of IcRn product, i.e., Josephson 
frequency, on the level characteristic for tunnel SIS junctions. Moreover, it should be noted 
that this property of the considered structure is unique. In the conventional SFS devices 
in order to reach the 7r-state it is necessary to realize either Lp > (2 — 3)£p or very large 
values of the Jbfs ^ 1 parameter at the SF interfaces. In both cases the IcRn product in 
the 7r-state is strongly reduced^ - — 

External magnetic field. In the parameter range when SIsFS junction is in (la) mode and 
far from the — 7r transition, current-phase relation has a sinusoidal form, I(<p) ~ Isis sin (p. 
To calculate the dependence of Ic from external magnetic field, H ext , we may use the 
standard Fraunhofer expression, 

sin(7T$/$ ) 



Ic(H, 



ext — J- CO 



7T$/$ C 



(1) 
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where 

$ = W \L eff H ext + L F H N(n t - H )\ (2) 

is magnetic flux inside of SIsFS junction, $0 is flux quantum, L e ff = 2X L + L s + L F + Lj, 
Xl is London penetration depth of S electrodes, Lj is thickness of I layer, N is the full 
number of clusters, n-^ = N^^/N are the concentrations of clusters in the F layer oriented 
parallel (N^) or antiparallel (iVj,) to the direction of H ext in the saturation region and Hq 
is an average magnetic field generated by a single magnetic cluster. In our simple model, 
Hq is assumed to be a constant, while and are functions of H ext . We assume further 
that the probability density p(H ext ) of a flip of a cluster in H ext is described by a Gaussian 
distribution 




H ext (Oe) 

FIG. 5. Experimental dependence^ of critical current Iq versus increasing (open circles) and 
decreasing (open squares) external magnetic field H ex t. Solid and dashed lines present the micro- 
scopic fitting of the data. Inset shows the theoretical and experimental magnetization loops versus 
external magnetic field H ex f 

p(H ext ) = (V2/V^S) exp [- (H ext - H*f /26 2 } , (3) 

where H* is the value of magnetic field, at which the flip of a cluster takes place. The 
expectation, H*, and the standard deviation, 5, in ([3]) are independent on H ext values. 
These parameters, as well as HqN product (i.e. saturation of the magnetization magnitude) 
in Eq.(J2]) can be found by fitting the magnetization curve 4nM(H ext ). 

We apply this model to explain the data experimentally observed in Ref.— in SIsFS 
structures having cross-section area 10 * 10 fim 2 and F layer thickness 14 nm. Figure 5 
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demonstrates the experimental dependencies of critical current Ic versus increasing (open 
circles) and decreasing (open squares) external magnetic field H ext . Solid and dashed lines 
present the microscopic fitting of the data. From hysteretic dependencies of F layer mag- 
netization shown in insert in Fig. we get H* « 11.4 Oe, 5 ~ 13 Oe and H Q N product 
~ 163 G. Remaining fitting parameters can be set from the current magnitude in the main 
maximum Icq ~ 2.9 mA and from the difference between zeros of I C {H) dependence in 
the saturation region, where oscillation period depends entirely on the effective length of 
structure L e // ~ 150 nm. 

The initial strong magnetic field H ext = — oo remagnetize all clusters of F layer into the 
homogeneous state n^(H ext ) = 1. Gradual growth of the external field provides the conven- 
tional Fraunhofer pattern (solid line) with expected maximum at the positive value H ext cor- 
responding to zero flux $ = 0. However, the clusters start to flip around the point H ext = H*, 
n^(H ext ) = J p(H ext )dH ext . As a result, the period of Fraunhofer oscillations decreases. 

— oo 

Similar situation takes place during field sweeping in the opposite direction (dashed line), 
from large positive to negative values. The densities n^(H ext ) and n^(H ext ) = 1 — n^(H ext ) 
can be described by the expression, n^(H ext ) = 0.5 (l + erf ((H ext =)= H*)/\/25)) , for for- 
ward and backward remagnitizations, respectively. Here erf(x) is the error function. 

In this paper we have demonstrated a number of unique properties of SIsFS Josephson 
junctions. These structures exhibit a large IcRn product in the 7r state comparable to that in 
SIS tunnel junctions commonly used in SFQ devices. Moreover, the whole structure behaves 
as a single junction with respect to an external magnetic field H ext . Based on that, we have 
developed simple model describing the behavior of the critical current in these junctions in 
external field H ext taking into account remagnetization of the F-layer. The model explains 
asymmetric Fraunhofer oscillations Ic{H ext ) in Nb-Al/AlO^-Nb-Pdo.ggFeo.oi-Nb junctions 
reported in^—. These effects provide the possibility to realize magnetic memory compatible 
with energy-efficient SFQ digital circuits^ with high switching speed. 

This work was supported by the Russian Foundation for Basic Research, Russian Ministry 
of Education and Science, Dynasty Foundation, Scholarship of the President of the Russian 
Federation and Dutch FOM. 
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